close all;clear all;clc;

sigma=0.2;
%равномерное распр-ие
y=rand(1,1e5)*(22-10)+16;
figure (1);
hist(y,120);

%нормальное распр-ие
y2=1*0+sigma*randn(1,1e5);
figure (2);
hist(y2,1200);

%хи-квадрат
k=2;%степень свободы
y3=sum(randn(k,1e5).^2);
figure (3);
hist(y3,1200);

%Релей
y4=sigma*sqrt(-2*log(1-rand(1,1e5)));
figure (4);
hist(y4,1200);

%Райс
x=0.38*1+sigma*randn(1,1e5);%+1;
y5=sqrt(x.^2+y2.^2);
figure (5);
hist(y5,1200);

%Холецкий
qcno_dB=45;
N=1e4;
stdn=8;

stdn_IQ=sqrt(stdn^2*N/2);

ro1=0.75;
ro2=0.5;
Dp=stdn_IQ^2; Dpe=ro1*stdn_IQ^2;Del=ro2*stdn_IQ^2;

L=chol([Dp Dpe Dpe;
        Dpe Dp Del;
        Dpe Del Dp])';
    
Nj = 1e6;
nIp=nan(1,Nj);
nIl=nan(1,Nj);
nIe=nan(1,Nj);
for j = 1:Nj
    nI = L*randn(3,1);
    nIp(j) = nI(1);
    nIe(j)=nI(2);
    nIl(j)=nI(3);
end





